B -> Xd(lP, %P)K decays in QCD factorization and X(3872) 
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B — * XciO-P>2P)K decays are studied in QCD factorization by treating charmonia as nonrela- 
tivistic bound states. No infrared divergences exist in the vertex corrections, while the logarithmic 
end-point singularities in the hard spectator corrections can be regularized by a momentum cutoff. 
Within certain uncertainties we find that the B — > Xci(2P)if decay rate can be comparable to 
B -> X ci(lP)K, and get Br(B° — » x'c\K°) = Br(B + — > x'ciK + )~ 2 x 1CT 4 . This might imply a pos- 



sible interpretation for the newly discovered X(3872) that this state has a dominant J PC 
cc component but mixed with a substantial D°D*° + D*°D° continuum component. 
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PACS numbers: 12.38.Bx; 13.25.Hw; 14.40.Gx 

The naively factorizable decay |l[ B — > \ c iK was 
studied|2| in the QCD factorization approach || in which 
the nonfactorizable vertex and spectator corrections were 
also estimated, but the numerical results were four times 
smaller than experimental data. Recently, these de- 
cays were also studied in the PQCD approach £|. In 
both the above approaches, light-cone distribution am- 
plitudes(LCDAs) were used to describe Xci- As argued 
in Ref. [5j, a more appropriate description of charmo- 
nium is the nonrelativistic (NR) wave functions which 
can be expanded in terms of the relative momentum q 
between charm and anticharm quarks. This argument 
is based on the nonrelativistic nature of heavy quarko- 
nium [Tj. With careful studies, we find that the two 
descriptions (i.e.LCDAs and NR) are equivalent for the 
S-wave charmonium states (see, e.g. [6|), but in the case 
of P-wave states the light-cone descriptions lose some im- 
portant contributions in the leading-twist approximation. 
This is not surprising since q can be neglected in S-wave 
states, but cannot be neglected for P-wave states even in 
leading order approximation. 

On the phenomenological hand, the study of B — > 
Xci{2P)K may help clarify the nature of the recently 
discovered resonance X(3872) |8j, since the measurements 
for A(3872) favor J PC = 1++ 9] and hence x cl (2P) be- 
comes one of the possible assignments for it. On the other 
hand, aside from the conventional charmonium | id ll l| . a 
loosely bound S-wave molecule of D°D*° + D*°D° has 
been suggested for X(3872)0,E1- 

Motivated by the above considerations, in this pa- 
per we study the decays B — > Xci(l-P 2P)K within 
the framework of QCD factorization by treating the 
charmonia % c i(lP, 2P) as nonrelativistic bound states 
with m c /mb taken to be a fixed value in the heavy 
b quark limit. We will estimate the production rate 
of Xci(2P) and argue that the X(3872) may be domi- 
nated by the \ci (2P) charmonium but mixed with some 
D°D*° + D*°D° continuum component. 

In the non-relativistic bound-state picture, charmo- 
nium can be described by the color-singlet NR wave func- 
tion. Let p be the total momentum of the charmonium 



and 2q be the relative momentum between c and c quarks, 
then v 2 ~ Aq 2 /p 2 ~ 0.25 can be treated as a small expan- 
sion parameter pj . For P-wave charmonium Xci j because 
the wave function at the origin 7?.p(0) = 0, which corre- 
sponds to the zeroth order in q, we must expand the 
amplitude to first order in q. Thus we have 

x^g)Tr[O"(0)P ls >,0)+O(0)Pr s >,0)], (1) 

where 0{q) represent the rest of the decay amplitudes 
and Pis^ (p, q) is the spin-triplet projection operator, and 
O a , P a stand for the derivatives of O, P with respect to 
the relative momentum q a The amplitudes O(q) can 
be further factorized as product of B — > K form factors 
and hard kernel or as the convolution of a hard kernel 
with light-cone wave functions of B meson and K meson, 
within QCD factorization approach. 

After q° is integrated out, the integral in Eq. is pro- 
portional to the derivative of the P-wave wave function 
at the origin by 



-^«m(«) = -ie* a (L z ),Pt- n' P (0), (2) 



(2tt) 3 



where e a (L z ) is the polarization vector of an angular 
momentum-1 system and the value of 1Z P (0) for char- 
monia can be found in, e.g., Ref. [l5|. 

In contrast to the NR description of Xci , the K-meson 
is described by LCDAs 0: 



l Jf f dxe^v'^+vp'-^U'^tKiy)} , (3) 
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where y and y = 1 — y are the momentum fractions 
of the s and d quarks inside the K meson respectively, 
and <Pk(x) is the leading twist LCDA of K-meson. The 
masses of light quarks and K meson are neglected in 
heavy quark limit. 
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regularization(NDR) scheme are given by 
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FIG. 1: Feynman diagrams for vertex and spectator correc- 
tions to B -^XcoK. 



The effective Hamiltonian for B —>XciK reads ^(| 



Hes^ (VcbV^CxOx + C 2 2 )- V tb V t * CiOi) , (4) 

i=3 



where Gf is the Fermi constant, Ci are the Wilson co- 
efficients and V qiq2 are the CKM matrix elements. The 
relevant 4-fermion operators Oi can be found in p|. 

According to all nonfactorizable corrections are due 
to Fig^ These corrections, with operators Oi inserted, 
contribute to the amplitude 0(q) in Eq. where the 
external lines of charm and anti-charm quarks have been 
truncated. Taking nonfactorizable corrections in Fig^ 
into account, the decay amplitude for B —>XciK i n QCD 
factorization is written as 



N c 4-7T N c V /j, 
N r 4tt N r \ Li 



«5=C5 + 



C 6 a s Cf 
N n 4tt N c 



C 6 (-6+121n^ + / J +/ JI ), (6) 



where C F = (N% - 1)/(2N C ) and ,u is the QCD renormal- 
ization scale. 

The function // is calculated from the four vertex cor- 
rection diagrams (a, b, c, d) in Fig^and reads 



fi = 
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We find that the infrared divergences are canceled be- 
tween diagrams (a) and (b), (c) and (d) respectively in 
Figffl On the other hand, this function is different from 
that in Eq. (11) of Ref. even when a nonrelativistic 
limit wave function cf)™(u) — S(u — 1/2) is adopted, as 
we have mentioned. 

For the two spectator correction diagrams (e, f) in 
Figffl the off-shellness of the gluon is natural to be associ- 
ated with a scale fih, ~ ^/to;,Aqcd, rather than fih ~ wib- 
Following Ref. |3|, we choose /i = Vm^A/, as 1.4 Gev with 
A/i = 0.5 Gev in calculating the hard spectator function 
/// and then, in the leading twist approximation, we get 



/// = 



tts{Vh)Ci{Hh) 87T 2 f K f B 1 
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(8) 



where £ is the momentum fraction of the spectator quark 
in the B meson and Cj(/Xfe) (i = 1,4,6) are the NLO 
Wilson coefficients which can be evaluated by the renor- 
malization group approach |T(| . 

The spectator contribution depends on the wave func- 
tion tfis through the integral 
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12^1 —TZ P (0)e* ■ PbF^M^), (5) 



where e is polarization vector of xa E3- Here F% is 
the B — > K form factor and we have used the rela- 
tion F (Af 2 )/Fi(M 2 ) =1-2 0, with z = M 2 /m 2 B m 
4m 2 /m 2 and M is the mass of Xci, to simplify the struc- 
ture of ©. 

The coefficients a, (i = 2, 3, 5) in the naive dimension 



o S 
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Since </>b(£) is appreciable only for £ of order Aqcd/tib, 
As is of order Aqcd- We will choose As ~ 300 MeV in 
the numerical calculations 

If we choose the asymptotic form of the K meson twist- 
2 LCDA , 4>k{v) — 6y(l — y), we can find logarithmic end- 
point singularities in Eq. (JH) just like that in Ref. Q, and 
we parameterize it in a simple way, 
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Xd(3511) 
Xci(3953) 
Xci(3872) 



0-2 



0.199-0.051i 0.000+0.002i 0.004-0.002i 
0.247-0.042i -0.002+0. OOli 0.007-0.002i 
0.236-0.044i -0.002+0.001i 0.006-0.002i 



TABLE I: The coefficients m of B -» XdO-P, 2P)K with dif- 
ferent choices of ■ 

X c i 



The mass of x&O-P) M Xcl — 3.511 Gev is known, 
but the mass of the missing charmonium \ c i(2P) has 
to be estimated by, say, potential models. We choose 
M xii =3.953 Gev following Ref. [H|. Then the form factor 
F\{M 2 ) can be determined by light-cone sum rules pflj . 



Fi(A^J = 0.80, Fi(M*, 



1.14. 



(11) 



We also choose M- 



3.872 Gev and FUM*, ) = 1.06 

to study if the AT(3872) behaves like a Xc i(2P) in their 
b-production processes. 

For numerical analysis, we use the following input pa- 
rameters : 



: 160 MeV, 
= V0TGeV 5/2 , 



m b = 4.8 GeV,m B = 5.28 GeVJ K = 
f B =216 MeV 21], ft' 1P (0) =ft 2P (0) 
Ci{n) = 1.21(1.082), C 2 (/x) = -0.40(-0.185), 
C 3 (m) = 0.03(0.014), C 4 (m) = -0.05(-0.035), 
C 5 (m) = 0.01(0.009), C 6 (/i) = -0.07(-0.041), 
a s (/i) = 0.35(0.22). (12) 

In 1121) the /i-dependent quantities at |Hh, = 1.4 Gev (/x = 
4.4 Gev) are shown without (with) parentheses. 

Using the above inputs, we get the results of coeffi- 
cients which are listed in Table. [I] With the help of 
these coefficients a*, we calculate the decay branching 
ratios of decays B — > X ci(lP, 2P)K with two different 
choices of M x < and get 



Br(B°-> Xcl (3511)iO 
Br(P°^ X ^(3953)iO 
Br(P°-> x ^(3872)*: ) 



1.79 x 10 , 
1.81 x 10~ 4 , 
1.78 x 10~ 4 . 



(13) 



Our prediction of Br(P° -> x c i(3511)i^ ) is about 2 
times larger than that in Q , although it is still about 
two times smaller than the recent data [23] • The differ- 
ence between the theoretical predictions and experimen- 
tal data may not be as serious as it looks like if we take 
into account the following uncertainties: (i) We have used 
a moderate value of 1Z 1P (0) predicted by different poten- 
tial models |l5j in our calculation, and a larger value of 
1Z 1P (0) may enhance our prediction in Eq. I|13l) signifi- 
cantly, (ii) In evaluation of ///, we only use the leading 
twist LCDAs of K-meson, and large uncertainties will 
arise from the chirally enhanced higher twist effects • 
(iii) Since the squared velocity v 2 of the charm quark 
in charmonium is about 0.25-0.30, the relativistic correc- 
tions may be important for these decays. 



Note that although the form factor in (|11|) and the 
coefficient ai in Table. Q] increase evidently as the char- 
monium mass increases, the decreased phase space and 
kinematic factors in JSJ will make a balance, and result in 
similar decay branching ratios in the charmonium mass 
region 3.51-3.95 GeV, as shown in (fT5|) . If we neglect the 
order a s corrections (i.e., in the naive factorization Q), 
the ratios among these three branching fractions in l|13|) 
would become 1 : 0.74 : 0.69. As estimated in l|13|) . the 
branching ratios for Xc i(2P) are 



Br(P°- 
Br(P+- 



w 2 x 10" 
= Br(P°- 



>Xci^°)- (14) 



Comparing Eq. I|14|) with the measured channel of the 
X(3872) 0: 

Br(B + -> XK+) x B x = (1.3 ± 0.3) x 10" 5 , (15) 

we see that the produced X(3872) looks like the Xc i(2P) 
if Bx is sufficient small, say, 3 ~ 7%. A similar conclusion 
has recently been obtained in a comprehensive analysis of 
X(3872) production at the Tevatron and B-factories [23l] . 
On the other hand, if X(3872]_ is a loosely bound S-wave 
molecule of D D*°/D*°D [Hill, a model calculation 
gives a smaller rate [l4j compared with Eq. I|14f> : 



Br(P+ -> XK+) = (0.07 - 1) x 10" 



(16) 



which requires a larger Bx > 10% to be consistent with 
the experimental data (|15[) . They also predict: 

Br(B°->X(3872)K°) < 0.1Br(B + ^X(3872)K + ).(17) 

So the measurement of B x and Br (B° ^X (3872) K°) is 
very helpful to identify the nature of X(3872). 

Recently, a preliminary result for a new decay mode 
X -> D°D°7r was found by Belief: 



Br(B -> XK) x Br(X 

= (2.2 ±0.7 ±0.4) x 10 



D°D 7T Q ) 
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(18) 



Eq. dJ implies that B x < 10%, if it can be confirmed by 
further measurements. This would disfavor the sugges- 
tion that the X(3872) is a loosely bound S-wave molecule 
of D°D* /D* D° with predictions of both decay[l3| and 
production^]. 

The above discussions about the X(3872) is based on 
the assumption that the X(3872) is a pure charmonium 
Xc i(2P) state. But this cannot be the case due to the 
coupled channel effects and X(3872) being in extremely 
close proximity to the D D*°/D*°D Q threshold. Per- 
haps a more realistic model for the X(3872) (for further 
discussions see |25j) is that the X(3872) has a domi- 
nant J PC = 1 ++ (2P) cc component which is mixed with 
a substantial real D D*°/D*°D continuum component 
(the D + D*~ I D*~ D + continuum component is kinemat- 
ically forbidden to be mixed in X(3872) and it is the u — d 
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quark mass difference that causes this isospin violation) . 
Thus X(3872) will have the following features. (1) The 
production of X(3872) in B meson decays is mainly due to 
the J PC = 1 ++ (2P) cc component, as discussed above. 
The production of X(3872) at the Tevatron is also due 
to this cc component and associated higher Fock states 
containing the color-octet cc pair and soft gluons. As 
was argued [ll] for the prompt charmonium production 
that cross sections of D-wave charmonia (which were sug- 
gested as a tentative candidates for X(3872) in 11]) could 
be as large as J/ip or ip(2S) due to the color-octet mech- 
anism, the P-wave (2P) charmonium could also have the 
comparable production rate to J/ip or ip(2S). But this 
does not seem to be obvious for a loosely bound S-wave 
molecule of D D*°/D*°D . (2) On the other hand, the 
D D*°/D*°D continuum component in X(3872) will be 
mainly in charge of the hadronic decays of X(3872) into 
D D*°/D*°D or D°D°tt as well as J/ipp° and J/ipcj. 
The latter two decay modes (J/ipp° and J/ipuj) may come 
from the first decay mode D°D* a /D*°D° and a subse- 
quent rescattering final state interaction and therefore 
have the same decay amplitudes [A(J/ipp°)=A(J/ipu})] 
that are smaller than the first decay mode amplitude. 
(3) A substantial D°D* Q /D* Q D Q continuum component 
in X(3872) may reduce the production rates in Eq. ((Hfl . 
and will also reduce the A(3872) — > Jhp r ) decay width, 
which can be as small as 11 KeV [ly (note that this 
2P-1S El transition is sensitive to the model details, see. 
e.g. P3). This is much smaller then the hadronic decay 
widths. But a large rate for Xci(2P) -> 7'0(2S , )=6O-1OO 
KeV will be expected. These qualitative features are use- 
ful in understanding the nature of X(3872) and should be 



further tested and studied experimentally and theoreti- 
cally. 

In summary, we study the decays B — > \ c i(lP, 2P)K 
in QCD factorization by treating charmonia as nonrela- 
tivistic bound states. We find that there are no infrared 
divergences in the vertex corrections, and the logarithmic 
end-point singularities from hard spectator interactions 
can be regularized by a momentum cutoff. Within cer- 
tain uncertainties we find the B — > Xci(^P)K decay rate 
can be comparable to B — > Xci(^P)K, and get Br(B a — > 
XciK ) = Br{B+ -> x'ci K+ ) ~ 2 x !0~ 4 - This might im- 
ply that the X(3872) has a dominant J_ PC = 1++(2P) cc 
component but mixed with some D°D*° + D*°D° con- 
tinuum component. The qualitative features of X(3872) 
are discussed and should be further tested and studied. 

Note. After this work appeared in h ep-p h/ 0506222 
we learned some new results from BaBar |26[: Br(B + ^ 
X(3872)*T+) < 3.2 x 10~ 4 ,i? = ^g^g|gg) = 
0.50 ± 0.30 ± 0.05 . We also note that a recent paper]!!] 
(hep-ph/0508258) obtained similar conclusions to ours 
for the X(3872). 
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